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1. Improved Computation of Fresnel 
Integrals 



In a previous paper [1], this author presented formulas 
for numerical computations of the Fresnel cosine and 
sine integrals, 
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to six significant figures by using the first three terms of 
the Taylor expansions ofC{x) and S{x) for Ixl ^ 0.6, the 
first three terms of the asymptotic expansions of the 
auxiliary functions /(x) and g(x) for Ixl > 3, and mod- 
ified rational approximations of fix) and gix) for the 



mid range. These formulas proved hard to use because 
they involve too many numerical constants. 

It was found subsequently that a simpler and more 
accurate method of computation can be based on a for- 
mula derived by Boersma [2], 
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where p„ and q„ are numerical constants tabulated in 
Ref. [2] and the notation has been changed in order to 
avoid confusion with symbols used elsewhere in this 
paper. On substituting t = 'UT^I2 and x = 'uu^l2 this is 
transformed into 
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, x>1.6, 
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/„ = (8/TTr"^^„, g„ = -(8/TTr"v»- 



(2c) 



Hence it follows by separation of real and imaginary 
parts and comparison with Eqs. (la) and (lb) that 



fix) = '2fnX- 



Kx) = 2 



gnX 



(2d) 



These give results with errors less than 6 X 10"'° for 
\x\ < 1.6 if the first 11 terms of the expansions are 
carried. 

Software and algorithms for computing Fresnel inte- 
grals in Fortran and C (not based on this paper) are also 
available on the Internet [4,5]. 



The required coefficients, as computed from Boersma's 
data and Eq. (2c), are listed in Table 1 . Equation (2d) was 
tested by computing sample values of C(x) and Six) and 
comparing them to the values tabulated in Ref. [3] to 
eight digits. The agreement was perfect, which is consis- 
tent with Boersma's statement that the error of Eq. (2a) 
is less than 5 X 10"'°. 

Table 1. Numerical values of f„ and g„ 



n 


/« 


& 





0.318309844 





1 


9.34626E-08 


0.101321519 


2 


-0.09676631 


-4.07292E-05 


3 


0.000606222 


-0.152068115 


4 


0.325539361 


-0.046292605 


5 


0.325206461 


1.622793598 


6 


-7.450551455 


-5.199186089 


7 


32.20380908 


7.477942354 


8 


-78.8035274 


-0.695291507 


9 


118.5343352 


-15.10996796 


10 


-102.4339798 


22.28401942 


11 


39.06207702 


-10.89968491 



2. References 

[1] K. D. Mielenz, J. Res. Natl. Inst. Stand. Technol. 102, 363 (1997). 

[2] J. Boersma, Math. Computation 14, 380 (1960). 

[3] M. Abramowitz and I. A. Stegun, eds., Handbook of Mathemati- 
cal Functions, (U.S. Gov. Printing Office, Washington, DC, 
1972). 

[4] D. W. Lozier and F. W. J. Giver, Numerical Evaluation of Special 
Functions in Mathematics of Computation, 1943-1993: A Half 
Century of Computational Mathematics, W. Gautschi, ed.. Am. 
Math. Soc. (1994) pp. 79-125. Also available at http:// 
math.nist.gov/nesf. 

[5] R. F. Boisvert, The Architecture of an Intelligent Virtual Mathe- 
matical Software Repository System, Mathematics and Comput- 
ers in Simulation, Vol. 36 (1994) pp. 269-279. Also available at 
http://math.nist.gov/gams. 

About the author: Klaus D. Mielenz is a physicist and 
retired Chief of the Radiometric Physics Division of the 
NIST Physics Laboratory. The National Institute of 
Standards and Technology is an agency of the Technol- 
ogy Administration, U.S. Department of Commerce. 



The above method is implicitly contained in Ref. [3], 
which mentioned Boersma's paper as well as the rela- 
tionship between the integrals in Eqs. (la), (lb), and 
(2a). Boersma also gave an approximation formula sim- 
ilar to Eq. (2a) for ImI < 4, or Ixl ^ 1.6. This was not 
used in this work because in this range it is simpler to 
compute C{x) and S(x) by using their Taylor expan- 
sions [3], 
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S{x) = ^s„ x'"*' 



So-- 
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